Wavefunction Collapse and Random Walk 



Brian CoUett and Philip Pearle 
Department of Physics, Hamilton College, Clinton, NY 13323 
(February 1, 2008) 

Wavefunction collapse models modify Schrodinger's equation so that it describes the rapid evo- 
lution of a superposition of macroscopically distinguishable states to one of them. This provides 
a phenomenological basis for a physical resolution to the so-called "measurement problem." Such 
models have experimentally testable differences from standard quantum theory. The most well de- 
veloped such model at present is the Continuous Spontaneous Localization (CSL) model in which 
a fluctuating classical field interacts with particles to cause collapse. One "side effect" of this in- 
teraction is that the field imparts energy to the particles: experimental evidence on this has led to 
restrictions on the parameters of the model, suggesting that the coupling of the classical field to the 
particles must be mass-proportional. Another "side effect," is that the field imparts momentum to 
particles, causing a small blob of matter to undergo random walk. Here we explore this in order 
to supply predictions which could be experimentally tested. We examine the translational diffusion 
of a sphere and a disc, and the rotational difi^usion of a disc, according to CSL. For example, we 
find that the rms distance an isolated lO^^cm radius sphere diffuses is « (its diameter, 5 cm) in 
(20 sec, a day), and that a disc of radius 2 ■ 10~^cm and thickness .5 • 10~^cm difi^uses through 
27rrad in about 70sec (this assumes the "standard" CSL parameter values). The comparable rms 
diffusions of standard quantum theory are smaller than these by a factor iO~^*^. It is shown that 
the CSL diffusion in air at STP is much reduced and, indeed, is swamped by the ordinary Brownian 
motion. It is also shown that the sphere's diffusion in a thermal radiation bath at room temperature 
is comparable to the CSL diffusion, but is utterly negligible at liquid He temperature. Thus, in 
order to observe CSL diffusion, the pressure and temperature must be low. At the low reported 
pressure of < 5 ■ 10"^''Torr, achieved at 4.2°K, the mean time between air molecule collisions with 
the (sphere, disc) is Ri(80, 45)min. This is ample time for observation of the putative CSL diffusion 
with the standard parameters and, it is pointed out, with any parameters in the range over which 
the theory may be considered viable. This encourages consideration of how such an experiment may 
actually be performed, and the paper closes with some thoughts on this subject. 
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I. INTRODUCTION 



Schrodinger was troubled by the collapse postulate associated with Bohr's "Copenhagen" version of quantum theory. 
This requires a superposition of macroscopically distinguishable states (an ill-defined concept), upon observation 
(another ill-defined concept), to be suddenly replaced by one of those states. In his famous "cat paradox" paper 0| 
Schrodinger wrote this "is the most interesting part of the entire theory," saying that it prevented one from ascribing 
reality to the wavefunction "because from the realism point of view observation is a natural process like any other 
and cannot per se bring about an interruption of the orderly flow of events." Dynamical wavefunction collapse models 
resolve Schrodinger's difflculty, allowing one to ascribe reality to the wavefunction (somewhat ironically) by altering 
Schrodinger's own equation, so that the collapse takes place in orderly and well-defined fashion. 

The Continuous Spontaneous Localization (CSL) model [HH], based upon previous models by Ghirardi, Rimini and 
Weber (GRW) |^ and one of the authors g is the most well-developed collapse model at present [HQ]. In it, to 
Schrodinger's equation is added a term which contains a randomly fluctating classical field ^(x, t) that interacts with 
particles, bringing about collapse. 

Although collapse is the desired and main effect, there are also "side effects." Because collapse narrows wavefunc- 
tions, particles gain energy from the field in this process [p|-p^. Experimental tests [pT|-p^ have resulted in restrictions 
on the range of permissable parameters for the model, sugggesting that the coupling between the field and particles 
(which determines the particle's collapse rate) is proportional to particle mass: thus, for a material object undergoing 
collapse, its nucleons are much more responsible for this behavior than are its electrons. In this paper we discuss 
another side effect: the random impulses particles get from the field results in random walk of objects. Indeed, in one 
of the earliest attempts at a dynamical collapse model, Karolyhazy jl^ discussed such behavior. Here we discuss it 
in the context of the CSL model, in order to see if the effect is measureable. 
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We first consider a sphere undergoing translational random walk. Section II summarizes the needed formalism 
associated with the usual Brownian motion in both air at temperature T and in a radiation bath at temperature T. 
Section III summarizes the results (of calculations given in the appendices, as are most of the detailed calculations 
in this paper) associated with CSL-induced random walk of the sphere. Section IV puts numerical values into these 
equations in three realms of air-sphere interaction: viscous, molecular and impact, in order of decreasing air molecule 
number density. It becomes clear that, in order to observe CSL diffusion, the air density must be low enough so that 
the mean time between air-sphere impacts is large compared to the time over which diffusion may be observed. 

In Sections V and VI we turn to discuss respectively translational and rotational diffusion of a disc. Rotational 
diffusion of a disc will be the subject of our experimental proposal, because a small translation distance (e.g., the 
disc radius) becomes, when it is a distance of rotation, equivalent to a large fraction of 27rrad and therefore more 
readily discernable. We consider a disc rather than a sphere because a perfect homogeneous sphere undergoes no 
CSL rotational diffusion since its rotated quantum states are identical. (An actual sphere's rotated states are slightly 
different so it does undergo a very small amount of collapse and rotational diffusion.) We find, for example, that a disc 
of radius 2- lO^'^cm and thickness .5-10~^cm diffuses through 27rrad in about 70sec: this assumes the "standard" values 
of the two parameters, proposed by GRW |^ for their model and taken over into CSL, which may be characterized as 
the time A"-'^ = lO^^sec it takes an isolated nucleon in a superposition of two localized states separated by a distance 
greater than a — lO^^cm to collapse to one of those states. A pressure of < 5 • lO^^^Torr is attainable and, at 
this pressure, we find the mean collision time between air molecules and the disc is about 45 minutes. These results 
are so encouraging that we consider, in section VII, the full range of (A, a) parameter values over which the theory 
may be considered viable (as well as the parameter proposal of Penrose Q based upon gravity). This indicates that 
experiments to observe diffusion of small objects could provide a definitive test of CSL and other collapse models. 
Therefore, in section VIII, we make a preliminary experimental proposal whose details we hope to examine in a future 
paper. 



II. BROWNIAN MOTION REVIEW 



A. Diffusion 

It is useful to review the usual Brownian motion formalism p^ . The Fokker-Planck equation for the probability 
density /9(x, v, t) for the position and velocity of the center of mass (CM) of a randomly walking sphere (radius R, 
mass M, density D) in a thermal bath at temperature T is 



dp 
dt 



E 



j dp 1 dv^p (3 d'^p 



(2.1) 



where (as will be seen) t characterizes the time to reach thermal equilibrium and (2/3t)^/^ is the equilibrium rms 
diffusion distance in time t. Eq. (2.1) can be used to calculate averages: /(t) = J J (ix(ivp(x, v, t)/(x, v). By 
multiplying Eq. (2.1) by and integrating by parts, and likewise for , we obtain dx^ /dt = , dv^ /dt = —v^/t, so 
tiJ (0) exp —t/r and 



(assuming a;^(0) = 0). Likewise, dx^ /dt — 2x^v\ dx^v^ /dt — — x^v^ /t, dv^ /dt — —2v^ /t + 2P/t , so we obtain 



(2.2a) 



{Axf = x3 - xi = 2/3r 



t/T- 



1-e 



-tjr 



(2.2b) 



[Axf T > — ^ 1 

^ ' t«T 3r2 ^ 4t 



(Axf 



t » T 



2pt. 



(2.2c) 



We particularly call attention to the ^ t^ diffusion for t << r. 
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One may readily express the variables [3 and r in terms of physical quantities. From Mdv^ /dt — —Mvi/r we see 

that the damping force is —{M/t)v = — fv. According to the equipartition theorem, the equilibrimn value of v^'^ is 
kT/M so, by Eq. (2.2a), kT/M = /3/r. Thus 

P = kT/t T^M/i. (2.3) 



B. Viscosity Factor 

We now consider various expressions for ^. In the case of a sphere in a fluid, as is well known, according to Stokes, 
the drag coefficient is 

e = 67r77i?, [l„, « i?), (2.4) 

where rj is the viscosity of the fluid and Im is the molecular mean free path. 

If the fluid is a gas, as the gas density decreases there are three realms for air-sphere interaction. At high enough 
density so that Im « R, is the viscous realm: here ^ is given by Stokes law (2.4). At lower density, where Im >> R, 
is the molecular realm. Here a colliding molecule may be considered to have a thermal velocity distribution since its 
last collision before hitting the sphere occurs so far away from the sphere that it is unaffected by the sphere's velocity 
(i.e., viscosity and hydrodynamic considerations are irrelevant). However, in this case, many molecular collisions occur 
over the shortest resolvable time interval so the Brownian motion assumptions still apply, and ^ is given by Eq. (2.5) 
below. At very low density, over time intervals where individual molecular collisions with the sphere can be resolved, 
is the impact realm where the Brownian motion assumptions no longer apply. 

In the molecular realm. Stokes law is no longer accurate. Experimental investigations into the correction to Stokes 
law, begun by Millikan and continued to this day, are fitted by 

^ Gnr/R 



l + (;™/i?)[a + /3exp-(7i?//,„)] 

where e.g., recent measurements jlj] on polystyrene spheres give a ~ I, P ~ .6 and 7 w 1. In the limit >> i? it is 
readily shown [p^ , po[ , assuming specular reflection of air molecules (other assumptions moderately alter the numerical 
coefficients), that a — 3/2, P — 0. In this case, using 77 = (l/i)nmgulm [n is the gas molecular number density, mg 
is the mass of a gas molecule and u is its mean velocity) in the above equation, the dependence upon Im disappears 
as one expects, resulting in 

^ = (47r/3)nmgui?2 = (8/3)ni?2(27rmgA:T)i/^ (Z,„ >> R) (2.5) 

where we have used u — (SfcT/Trmg)^/^. 

In the case where the sphere moves in thermal radiation it is shown in Appendix D that a result of Einstein and 
Hopf ||2^,2^ may be adapted to obtain, for a dielectric sphere of large dielectric constant (but also true up to a 
numerical constant for other shaped objects, where R^ is replaced by the object's volume). 

When Brownian motion assumptions apply, regardless of the physical source of ^, the rms diffusion distance Ax's 
long time and short time behaviors differ. Einstein's well known result for t » t and the result for t << r follow 
from Eqs. (2.2c), (2.3): 

\2kTtY'^ 

Ax —- (2.7a) 



Ax 



t « T 



2kT^t^ 
MP 



1/2 



(2.7b) 
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III. CSL RANDOM WALK OF A SPHERE 



A. Diffusion of Center of Mass 



In the case of CSL, we consider an ensemble of sphere CM wavcfunctions (q|V', t)w, each evolving under a particular 
sample field «;(x, t). They are described by the density matrix p{t) whose evolution equation (see Appendix A) 
satisfies 



d{^\p{t)\<i') 
dt 



= -m 



.pit) 



2M 



|q') 



y2 



li 



dzdz' 



$(z - z') - $(z - z' + q - q') 



(q|pW|q') 



(3.1) 



under the approximation that the mass in the sphere is uniformly spread throughout it. In Eq. (3.1), P is the CM 
momentum operator, V is the volume of the sphere and $(z) = exp— z^/4a^. The electrons have been neglected 
because of their smaller mass and lower collapse rate, and the proton and neutron masses are taken to be equal for 
simpUcity, so Eq. (3.1) depends just upon the nucleon number N. 

Eq. (3.1) may be used to calculate ensemble averages of expectation values of operators: {F)(t) = Tr[Fp(t)]. The 
trace of the CSL term in Eq. (3.1) multiplied by the CM position operators (or any function of them), by 
or by Q^P^ + P^Q^ vanishes. Thus we obtain d{Q3)/dt = JP^/M, d{Pi)/dt = and so (Q^(i) = (assuming 
(QJ)(0) = 0) and JP^it) = (assuming (P^)(0) = 0). 

However, the trace of the CSL term in Eq. (3.1) multiplied by P^ docs not vanish. Collapse increases energy 
because it narrows wavefunctions and the references given in section 1 show that (neglecting the collapse behavior 
associated with the electrons) the rate of energy increase is given by 



3Ar7V2 



4Ma2 ' 

As is shown in Appendix A, for the CM part of the energy it follows from Eq. (3.1) that 



(3.2) 



d (P^^) Xh'N^f{R/a) 



dt 2M 



AMa? 



(3.3) 



The factor / essentially characterizes the collapse rate when the sphere is displaced by the distance a (see the discussion 
after Eq. (AlO) in Appendix A). f{R/a), given in analytic form in Eq. (A9b), is a monotonically decreasing function 
of its argument, with /(O) = 1, /(I) = .62 and f{R/a) 6(a/i?)"^ for R » a. Summing Eq. (3.3) over the three 
values of j and comparison with Eq. (3.2) shows that, for small R/a, the excitation of the CM accounts for almost 
all of the sphere's energy increase but, as R/a increases, internal nuclear excitation (too small to observe at present) 

accounts for more of it. 

Therefore, using Eq. (3.1), since d{Q^^)/dt = {Qi Pi + PiQi) /M , d{Qjpj + PiQi)/dt = 2{Pj^)/M and d{Pj'^)/dt 
is given by Eq. (3.3), we find 



9 I PH, 



>(o) 



Xh^f{R/a)t^ 



(3.4) 



where m is the mass of a nucleon. We note the ~ t^ diffusion associated with a random force without damping. 
This occurs essentially because the average square velocity is increasing so the distance of each "step" in the random 
walk increases with time. In Eq. (3.4) we have utilized M = Nm to emphasize that, for R « a, the diffusion is 
"universal," i.e., independent of the material and size (or, it turns out, shape) of the piece of matter and, in general, 
that the dependence on N is only indirect, through the sphere's radius R. 



B. Wavepacket Width 

Eq. (3.4) is the result needed to describe CSL random walk. However, it is necessary to show that the ^ t"^ term 
is not due to an increase in the width of the CM wavepackets in the ensemble _but truly due to th e diffusion of the 
centers of the packets. That is, the ensemble mean square wavepacket width is = {[Q^ — {Q^)]'^), what we want is 
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the ensemble mean of the squared displacement of the center of the wavepackets (Q-')^ and what we've got from Eq. 
(3.4) is {Q3'')=l^ + ^Qjf. 

Under the combined influence of the collapse (which tends to narrow wavefunctions) and the normal Schrodinger 
evolution of a free object (which tends to expand wavefunctions), tends to an equilibrium size in a characteristic 
time Ts- This has been discussed in the context of the GRW model [p| j2^ , |25t and for a particle in a simple continuous 
collapse model by Diosi 1^]. We discuss it for the CSL model in Appendix B. It requires a separate treatment because 

(Q^)^ (which is not (Q^) = 0) and so cannot be found from the density matrix since they involve an ensemble 
average over a quantity that is quartic in the statevector. 

According to Appendix B, the asymptotic CM wavepacket width is the same for every wavepacket (i.e., no ensemble 
average need be involved): 



t » T, 



2XmN^f{a/R) 



1/2 



(3.5) 



This expression for Soo can be understood as follows. If a wavepacket has width s, due to its Schrodinger evolution 
it expands a distance ~ {ti/Ms)At in time At. Due to the collapse evolution it contracts a distance ~ (collapse 
rate) (fractional decrease)s. As discussed in Appendix A (after Eq. (10)) the collapse rate is XN^f. The fractional 
decrease is that which would occur if a gaussian of width s << a is multiplied by a gaussian of width a, namely 
{s/ay^ . Thus it contracts a distance XN'^f{s/a)'^sAt. Equating the Schrodinger expansion to the collapse contraction 
and solving for gives the result ^(3.5). 
The characteristic time to reach this width is 



Nms%, 



(3.6) 



This can be understood as the time it takes a packet of width Sqo to expand by the distance s^a-. {fi/Msoa)Ts — Soo- 

Appendix B also shows that the CSL diffusive behavior soon becomes the dominant contribution to {Q^"^) once 
equilibrium has been reached: with initial equilibrium conditions, Eq. (3.4) becomes 



2r2 



t' 



12r3 



(3.7) 



Eqs. (3.5), (3.6) are derived in Appendix B under the assumption that s^{t) « a^. From Eq. (3.5) this implies 
> 3 - 10'' nucleons which, for ordinary matter densities (Igm/cc < D < 20gm/cc) means that Eqs. (3.5), (3.6) hold 
for i? > 2 • 10-^cm. 

IV. TRANSLATIONAL DIFFUSION OF A SPHERE: NUMERICAL VALUES 

We shall now put numbers into these equations so as to consider the conditions necessary to observe CSL-induced 
diffusion of a sphere. In the following we shall use the GRW parameter values, = lO^^sec and a = 10~^cm, until 
section VII when we consider the full range of allowable parameter values for the theory. 



IS 



A. CSL Diffusion Alone 

According to Eq. (3.4), under the CSL mechanism acting alone, the rms distance along an axis the sphere diffuses 

1/2 



AQJ = — 

ma 



6 



6.5/^/2 (i days)^'^ cm 



3/2 



(4.1) 



where / = 1 for R « a. Since Eq. (4.1) does not depend directly upon N it is independent of the density D. It 
does depend upon R through /, decreasing rapidly as R increases: 



AQ3 



R » a 



l%{a/Rf{t days)^^^ cm 



Table 1 lists AQ for various values of R and t. 
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Although diffusion is all wc shall be concerned with in the remainder of this paper, we give here some values for 
Soo and Tg. From Eqs. (3.5), (3.6) it follows that 



'R«a 



3.8- 10-^D-3/4(a/i?)9/4cm, 



R» a 



2.4- 10-^£)-3/4(a/iJ)5/4cm 



' R«a 



.58D-i/2(a/i?)3/2sec, 



' R»a 



.23D-'/^{R/ay/^sec 



4 • 10 ''cm in 



{D is in gm/cc). Tabic 2 lists values of s^c and for various values of R, for a sphere of density D = Igm/cc. (Wc 
note that Tg increases as R moves away from a in either direction since as R decreases the collapse rate decreases and 
as R increases the Schrodinger spreading rate decreases). 

For example, an i? = 10~^cm sphere's center of mass wavefunction reaches equilibrium size Sot 
Tg « .6sec, and diffuses AQ « (60microns, 5cm) in (lOOOsec, Iday). 

It is worth examining the diffusion to be expected were the CSL hypothesis to be false (A = 0) and the Copenhagen 
concept of collapse somehow occurring "upon observation" to be employed. In this case we utilize Eq. (3.4), AQqm = 
{P/M){0)t. If the sphere is observed at f = 0, localized to w 2R, we may take {P){0) ~ fi/AR. If the sphere is then "in 
the dark" (unobserved) until time t, wc obtain from AQqm ~ th/[D{4:/3)nR^AR], for R = 10"^cm and D = Igm/cc, 
that AQqm « (6 • 10~^cm, 5 • 10~'*cm) in (lOOOsec, Iday). These numbers are smaller than their CSL counterparts 
by the respective factors (10~^, 10"''). 



B. Diffusion in Air 



The CSL diffusion distances in vacuum given in Table 1 are much reduced by the — damping due to collisions 
with air molecules. 

[Another effect of these collisions is that, as the air molecules collide with the sphere they become entangled with its 
states, increasing the effective collapse rate. This effect is complicated, depending upon the sphere's quantum state's 
differences of air molecule density in sized volumes surrounding the sphere. Because the air molecule density is 
much less than the sphere density and because the quantum states of the sphere which compete in the collapse "game" 
are so spatially close, I shall ignore this effect.] 

The Fokker-Planck equation for the combined CSL and Brownian diffusion in air, which replaces Eq. (2.1), is 



dp 

di 



3 

E 



dp ^ dv^p 



+ 



+ 



d^p 



The long time and short time diffusion expressions which replace Eqs. (2.7a,b) are 



{Axf 



t»T 



2kT 



(4.2) 



(4.3a) 



{Axf 



t«T 



2kT£, 
3M2 Qm^a? 



(4.3b) 



1. Viscous Realm 



First consider the viscous realm. At room temperature Tq and atmospheric pressure po, the mean free path of air (N2 
or O2) is Im ~ .6- 10~^cm. For ^ we use Stokes' law (2.4) or the corrected equation following it (needed for R = 10~^cm 
since then Im ~ R, which amounts to a 40% decrease in ^ if a = 3/2 and /3 = 0): with 77 w 2 • 10~^gm/cm-sec we have 

T = M/^ w (2 • 10-^ 10^)sec for R = (10-^ l)cm 



Then we may apply Eq. (4.3a) for t > r to obtain 



Axbr = 



2kT 



1/2 



(.6, 1.4 • 10"^) (Mays) ^/^cm for R = (10-^ l)cm 



(4.4a) 
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AxcsL = f y j ( 2^42 7 ^ ' 10""^(*days)i/2(,j^ foj. ^ > iQ-Scni (4.4b) 

(note that Eq. (4.4b) is independent of R). For t < t, where Eq. (4.3b) apphes, the ratio Axcsl/ Axbb. is the same 
as that given in Eqs. (4.4). 

Clearly the CSL diffusion is swamped by the Brownian diffusion in the viscous realm (especially since it is the sum 
of the squares of Eqs. (4.4) which add in Eq. (4.3a)). 



2. Molecular Realm 



We next turn to the molecular realm where Im » R, which can be achieved by lowering the air density through 
lowering the air pressure p (which we shall give in units of picoTorr: lpT=10~^^T). We focus upon R = 10~^cm 
spheres since, ii R >> a, Axcsl decreases ~ R~^ and so is less easily observed (and also Axbr ^ R^^ so no relative 
advantage is gained by increasing R). Moreover, spheres of this size have been typical of observations of Brownian 
motion in air ||l7| , p7|| . From Eq. (2.5) we find the time to reach thermal equilibrium is 

T = M/^ « 2 • 10^(r/ro)i/2(ppT)-isec 

Since r is so long for picoTorr pressure or less, we apply Eq. (4.3b) for t << r: 

1/2 

w 2 • 10-'^{ppTy^^iT/To)^/^D-H^/'^cm (4.5a) 
1/2 

w 2 • IQ-'^i^/^cm (4.5b) 
{t is in sec). 

It follows from Eqs. (4.5a,b), even at liquid He temperature T — 4.2°K and for a dense sphere D — lOgm/cc, that 
Axbr. ~ AxcsL requires the very low pressure of p « 10~'^pT. But, at this low pressure, the Brownian assumption 
of many molecule-sphere collisions occurring in the shortest observable time interval no longer applies. Therefore, 
observation of CSL diffusion requires the air density to be low enough to be in the impact realm. 



Axbr = 



2fcTe 3 
3M2 



AxcsL 



3. Impact Realm 

Consider the mean time between molecule-sphere collisions. The mean number of collisions/sec-area of air molecules 
is (nu/4) so the mean time between molecule-sphere collisions is Tc = [{riu/ A){AttR^)]~^ . Moreover, the change in 
speed of the sphere due to one collision with an air molecule is Av ~ u{mg/M). With u « 4.5 • 10''^(T/To)^/^cm/sec 
and n « 2.5 • 10^^{p/po){To/T)cm^^ wc obtain, for an i? = lO^^cm sphere, 

Tc « 2(T/ro)i/2(ppT)-isec, Av ^ 5 ■ 10-^{T/Toy/'^cm/sec. (4.6) 

[Incidentally, we can understand the molecular realm's Brownian motion in terms of the impact realm motion if we 
consider the impact realm but for t >> Tc, so that many collisions have occurred in time t and so Brownian motion 
considerations apply. Then Axbr in Eq. (4.5a) may be written in terms of the quantities in Eq. (4.6), using ^ taken 
from Eq. (2.5): 

Axbr ~ [ikT/M^)nR^{mgkT)^/HY/^ ~ Avt[t/T,]^/^ , 

This says that the Brownian diffusion distance in time t is the distance the sphere goes with the speed it gets from a 
single collision multiplied by the square root of the number of collisions (the expected fluctuation in the number of 
collisions)]. 

Eq. (4.5b) (the same as Eq. (4.1) gives the CSL diffusion distance in the impact realm. In order to observe CSL 
diffusion over the largest distance, one wants Tc to be as long as possible and so. By Eq. (4.6), one wants the lowest 
possible pressure. An experiment conducted at pressure< 5 • 10~^^Torr at 4.2°K has been reported In these 

conditions, the mean coUsion time is ~ 80min. In this time, according to Eq. (4.5b), Axqsl ~ .7mm. This should 
be readily observable, so much that it encourages one to contemplate an experiment to test CSL over a wide range of 
parameter values (Section VIII). 
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C. Diffusion in a Thermal Radiation Bath 



For completeness, we note that, even when one eliminates collisions of the sphere with air molecules over some 
sufhciently long time interval, there is still thermal radiation to supply damping and random impacts and thus induce 
Brownian motion. However, as we shall soon see, this is very small at liquid Helium temperature. 

The viscosity coefficient in the case of thermal radiation, obtained in Appendix D and cited in Eq. (2.6), has the 
numerical value 

^RAD - 4 • 10-29 (i?/10-^)6(T/To)8gm/sec (4.7) 

{R is in cm). 

It follows from this that the time (2.3) to reach thermal equilibrium is 

TRAD = M/^RAD ~ 10'^DiR/lO-YHT/To)-^sec. (4.8) 

At room temperature or less, trad is so long that only the case of t << trad is of interest. Then, using Eqs. (2.7b) 
and (2.6), 

1/2 

w &D-\T/Tof^^{t/10'^f/^cm. (4.9) 

(note that (4.9) is independent of R). According to Eq. (4.9), at room temperature a. D — Igm/cc sphere of any 
radius will diffuse « 7cm/day due to thermal radiation alone. This is essentially equal in magnitude to the CSL 
diffusion for an i? = lO^^cm sphere (which is unaffected by such a small damping coefficient). However, at T = 4.2°K 
we have Axrad ~ 4 • IQ-^cm in a day which is utterly negligible. 

Therefore, at liquid He temperature, which is needed to obtain the low pressure of the impact realm, we need not 
consider the random walk due to thermal radiation. 



V. TRANSLATIONAL DIFFUSION OF A DISC 

Because rotation through 27rrads may be more easily detected than a comparable translation, our experimental 
proposal (section VHI) is based upon observing rotational diffusion. However, since a uniform sphere displays no 
CSL rotational diffusion we consider a more asymmetrical object, a disc. In this section, for completeness, we discuss 
translational diffusion of a disc. In section VI we shall discuss rotational diffusion of a disc. 



A. Brownian Diffusion 

Consider a disc of radius L and thickness b. For Brownian motion, the time dependence of the rms diffusion is given 
in section II in terms of ^ (e.g., Eqs. (2.7)). In the viscous realm, for an oblate spheroid (cc^ +y'^)/L'^ + 2^/(6/2)^ = 1 
(close enough to a disc). Lamb ||2^ shows that, for b « L, 

^ « I67^L (motion perpendicular to face) (5.1a) 

^ « (32/3)77^ (motion along edge). (5.1b) 

This is not qualitatively different from Stokes law (2.4) (with Rk, L). 

In the molecular realm where Im » iL,b), one may readily calculate, as in references [ p^ , pO| , assuming specular 
reflection of the molecules in the disc rest frame, 

^ — 4:nL^{2TTmgkT)^^'^ (motion perpendicular to face) (5.2a) 

^ ^ 2nLb{2TmigkTf/'^ (motion along edge). (5.2b) 

Eq. (5.2a) is not qualitatively different from Eq. (2.5) for a sphere (with L k, R). As for Eq. (5.2b), decreasing b to 
reduce ^ for edgewise motion does not reduce Axbr/ ^^xcsl since, from Eqs. (4.3), [Axbr/ Axcsl]'^ ^ ^/A'P ^ b^^. 
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B. CSL Diffusion 



For CSL, the time dependence of the rms diffusion is given in Appendix A, Eq. (A. 10) (and copied in Eq. (3.4)). 
/ = 1 for a disc with aU dimensions << a, so in this case the disc's diffusion is no different from that of a sphere with 
R « a, Eq. (4.1). There is a difference for (6/2a)2 << 1 and {L/2af » 1: 

/ {2a/ L)^ (motion perpendicular to face) (5.3a) 

/ ^ (4/V7r)(a/L)^ (motion along edge). (5.3b) 

Thus the (thin) disc diffusion decreases less with increasing size than does the sphere's diffusion, for which / — > 
6{a/R)'^. Thus, if a larger object is needed for greater visibility, a larger radius disc gives greater diffusion than does 
a sphere of the same radius. However, the conclusion reached in section IV for a sphere holds as well for a disc: the 
impact realm is required to effectively remove Brownian motion in order to see CSL translational diffusion of a disc. 



VI. ROTATIONAL DIFFUSION OF A DISC 



Rotational Brownian motion was (naturally) first considered by Einstein The Fokker-Planck equation for an 
object rotating about a fixed axis through angle 9 with angular velocity uj is identical in form to Eq. (2.1) with the 



replacements 



For rotation, the viscous torque on a sphere is —£,r,ot^ where 
S.ROT = 87r77i?3, (/„ « i?). 



which replaces Stokes' law, Eq. (2.4). Eq. (2.7a) is replaced by 



BR 



t » Trot 



2kTt 



S,ROT 



1/2 



(6.1) 



(6.2) 



where trot = I/£,rot and / = (2/5) AIR"^ is the sphere's moment of inertia. We see from Eqs. (6.1), (6.2) compared 
with Eqs. (2.4), (2.7a) that AObr « Axbr/R- 

But, the case of a sphere is of no use to us. In the approximation we make, where the nuclear mass is uniformly 
spread out over the sphere, there is no difference between two rotated quantum states of the sphere so, according to 
CSL, there is no collapse and therefore no random rotational motion (without this approximation there is collapse 
and random rotation but it is very slow). However, CSL random rotation does occur for a nonspherical object. 



A. Brownian Rotational Diffusion 



Here we consider rotational diffusion of a disc "on edge" (i.e., oriented with the face of the disc in a vertical plane), 
of radius L and thickness 6 << L, in the molecular and impact realms (Im » L). For a sphere in these realms, if 
the molecules make elastic collisions with the sphere, they do not transfer momentum parallel to the sphere face and 
so do not cause any torque {^rot = 0). However, for the disc, a straightforward calculation (as in [|l9|,^) yields the 
torque = —^rot'-^ about an axis passing through the edge and center, where 

UoT = (4/7r)nL4(2^TO<,fcT)i/2 (6.3) 

and Eq. (2.7b) is replaced by 

'2kT^ROTt^ 



A9 



BR 



t << Trot 



3/2 



(pprmr/j^l^ 

~ 80— ; — , , , , — — Trrads (6.4 

(L>gm/cc)(M/i)(id^)2 ^ ' 



{I w ML'^/A). Here we have employed the rather weird unit Id/i — lO^^cm because the dimensions of the disc we 
are considering are such that the factors (bd^), (id/i) are of the order of unity. We only give Eq. (6.4), valid for 
t << Trot, because the time to reach thermal equilibrium is so long at picoTorr pressures or less: trot = 1/S.rot ~ 
5 • 10^15(6 in d/i)(ppT)-i(T/To)i/2sec. 

According to Eqs. (6.4) and (6.5) (below), Brownian motion dominates CSL diffusion even at IpT pressure, for 
discs with dimensions of the order of a = Id/i. Therefore one must go to lower pressure, to the impact realm, to see 
CSL rotational diffusion. 
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B. CSL Rotational Diffusion 



For CSL rotational diffusion, it is shown in Appendix C that 

" (%^) " -OlS/^oV/^-d (6.5) 

where FIG. 1 contains a graph of fROT{c(, vs. a = {L/2a) for various vahics of /3 = (6/2a). For example, 
/rot ~ 1/3 for b « .5a and L w 2a. 

For this example, according to Eq. (6.5), A9csl diffuses through 27rrad in about 70sec. If A were 10"'' times 
smaller i.e., A « 10~^° (and still a = 10~^cm), which is at the edge of where one may consider the theory to be viable 
(section VII, Eq. (7.3)), this time is about 25 min. 

It is worth examining the rotational diffusion to be expected from standard quantum theory (as was done for 
translational diffusion at the end of section IVA). From Eq. (C6) with A = we have 

aoqm = {cmt/i 

where (£) (0) is the expectation value of the angular momentum operator in the initial state. If the disc is observed at 
t = 0, localized to k. -k/A, then we may take (>C)(0) « ?i/[2(7r/4)] = 2h/iT. If the disc is "in the dark" (unobserved) 
until time t, using I = (Z)67rL^)(L^/4), we obtain 

With the choices D = Igm/cc, 6 = .5 • lO'^cm, L = 2 ■ 10-''^cm, we get AOqm ~ lO^-'^trad. Thus, A9qm »(.l, 
1, 86)rad in t =(100sec, lOOOsec, Iday). These numbers are smaller than their CSL counterparts by the respective 
factors (100, 300, 3000). However, were A sufficiently small, this diffusion could be observed in the experiment we 
propose (section VIII). 



C. Gas-Disc Collisions 



The times given in the example of the previous section for diffusion through 27rrad (« 70sec for A~^ = lO^^sec, 
f«25min for A~^ = 10^°sec) should be compared with the mean time between collisions of air molecules with the disc. 

Assume Nitrogen molecular gas at temperature 4.2°K and pressure 5 • 10~^^Torr. The mean molecular speed is 
u — [8A:T/7rmg]^/^ « 5.6 • lO^cm/sec. The molecular density is p = p/kT w 115particles/cc. The molecular flux is 
J = pu/4 « 1.5 • lO^particles/cm^-sec. 

The mean time between collisions (we consider collisions with the 2 faces of the disc but neglect collisions with the 
edge) is thus Tc ~ 1/(2 JttL^) w 45min. 

We conclude this subsection with an estimate of the effect of a collision. A Nitrogen molecule with speed u impacting 
perpendicular to the disc face at distance L from the rotation axis ("worst possible case") conveys to the disc an 
angular velocity 

m„uL 33 , , ,^ „, 

For the example we have been considering, this is w 8rad/sec. Such a sudden jump in the angular velocity should 
be readily observable and distinguishable from the expected CSL behavior. 



VII. PARAMETER VALUES 



In the previous sections of this paper, for clarity's sake, in numerical calculations we have used the values of the 
CSL parameters (A^^,a) suggested by GRW, namely (lO^^sec, lO^^cm). However, these values have no theoretical 
underpinning and were simply chosen to give reasonable results: other values are possible. However, not all values 
are possible. 

Therefore we examine already existing experimental and theoretical constraints on these parameters. Any new 
experiment must be considered as placing further constraints. While one may hope that experiment reveals an 
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"anomalous" random walk confirming the existence of a CSL-type collapse process and disclosing the values of the 
parameters, one should consider the possibility that this does not happen. We consider the additional constraints 
negative results could provide. In particular, we consider what would be needed to eliminate CSL as a viable resolution 
of the "measurement problem" . 

We also discuss two other topics. One is the random walk associated with a suggestion by Penrose of a connection 
between gravity and collapse. This can essentially be interpreted as giving the results of this paper with a particular 
value for A. We also show the range of parameter values consistent with a speculation that the fluctuating field w 
has a thermal basis, based upon cosmological considerations and an analogy between standard random walk and CSL 
random walk. 

In what follows, means up to a numerical factor not too far from 1. 

An experiment which looks for photons emitted by the atoms in an underground shielded slug of Germanium places 
a limit on the number of bound electrons or nucleons "spontaneously" excited in Ge atoms [|2|,^. Spontaneous 
excitation of bound states is expected from the CSL collapse mechanism, which narrows electron and nucleon wave- 
functions thereby giving these particles increased energy (presumably the energy for this comes from the fluctuating 
collapse-causing field w [0). For example, a Is electron ejected from an Ge atom will result in radiation of an 11.1 
keV (equal to its binding energy) shower of photons from the atom's remaining electrons as they cascade downward 
added to radiation equal to the kinetic energy of the ousted electron which it rapidly loses in collisions with other 
atoms. 

The present experimental upper limit on the rate of photon pulses appearing in 1 KeV bins above 11 keV is ~ 
.05 pulses/(keV kg day) [|2|. The theoretical excitation rate is conveniently expanded in a power series in (size of 
bound state/a)^. The first term in this series turns out to vanish identically if the collapse coupling constant is mass 
proportional ||ll|. We have assumed this in the present paper (e.g., see Eq. (A.l) et. seq.) because, for atomic 
spontaneous excitation, the numerical coefficient of this first term is large enough to make the experiment sensitive 
to the relative coupling constant size of electrons and nucleons, and the results make mass-proportionality likely. The 
experiment is less sensitive to the second term in the series but the data on excitation rate of nucleons still provides 
a constraint (because of the now-assumed small electron coupling constant, the electron excitation rate data does not 
provide as strong a constraint). The theory gives probability/sec^ A(nucleon diameter/a)^. Since, in Ge, there are 
8.3 • lO^^atoms/kg and A 72 and, using nuclear radiusw 1.4 • 10~^^A^/"^cm and 8.6 • lO^sec/day, we get 

A-ifl^ > 2-10-15. 

However, the strongest present constraint, based upon the same experimental data, was provided by Fu |]33|| . He 
calculated the rate of radiation by a free electron (mass me) due to being shaken by the collapse mechanism. He 
obtained for the number of photons of energy E radiated by an electron per second per energy the expression 

^ \{m Imfe^n ^ {X/a') f j^^y) _^ 2.1 • 10-«counts/(keV kg day). 

The last term on the righthand side of the above equation gives the rate of radiation from the 4 valence electrons 
(essentially free) from each atom in a slug of Gc (using 8.29 • lO^" atoms/kg for Ge) at « 11 keV with the GRW 
parameter values. This and the experimental upper limit quoted above leads to the experimental constraint: 

X-^a^ > .4 (7.1) 

This constraint (labelled line 1) is graphed in Fig. 1: the allowed region is to the right of the line. 

Diffusion experiments should do much better than (7.1) in constraining the parameter values. For example, consider 
a rotational diffusion experiment such as we sketch in the next section. One expects to be able to detect a AO « 7r/2 
diffusion in 45 minutes. If such a diffusion were not detected, Eq. (6.5) (where AQcsl goes as A^/^/a^) gives 

A-ifl-* > 10^ (7.2) 

(region to the right of the line labelled 2 in Fig.l). This amounts to being able to detect AQcsl a factor 10^^ times 
smaller than that expected using the GRW parameter values. 

There is also what may be called a "theoretical constraint" ||l^,|l^, although it is fairly rough. The purpose of a 
collapse model is to account for the world as we see it. The model may be considered to fail if it allows an observable 
object to remain in a superposition of two well-separated locations "too long". How long is "too long"? We might 
take that to be human perception time ~ .Isec. 

For a first example, consider an object which is just visible, a sphere of diameter 4 • lO^^cm, in a superposition 
involving a displacement >> a, with a > 4 • lO^^cm. The collapse time is ~ X^^/N"^, where N is the number of 
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particles in the sphere. If the sphere's density is D w Igm/cc, then N k, 2 ■ 10^° and the condition A ^/A^^ < .Isec 
imphes 

A-i<4-10i^ (7.3) 

(region below the line labelled 3 in Fig.l). 

For a second example, again consider the above sphere but with a superposition involving a displacement < a. In 
this case the collapse time is ~ (4A~"'^a^)/[A^-displacement]^. Using the smallest possible discernible displacement, 
4 • 10~^cm, the condition that the collapse time is < .Isec implies 

A-^a^ < 1.6- 10^° (7.4) 

(region to the left of the line labelled 4 in Fig.l). 

These "theoretical constrains" are rough but we may take them seriously enough to observe, from Fig. 1, that 
constraints (7.2) and (7.4) still permit a narrow wedge-shaped range of allowed parameters. 

However, suppose one were able to perform an experimental test of translational diffusion of a sphere with a precision 
for AQ that is 10~^ times smaller than AQcsl with the GRW parameters, and find a null result. According to Eq. 
(4.1) the resulting constraint is X^/^/a < 10-3(10"i^sec-i)i/2/(io-5cin) or 

X-'a^ > IQi" (7.5) 

It appears that the conflict between (7.5) and (7.4) would make CSL nonviable. 
We close this section with two additional considerations. 

First, in Appendix E we argue that a proposal by Penrose |6|, and other suggestions involving a gravitational basis 
for collapse ||3^-|3^, arrive at an effective value for A: Xg ~ Gm? /ah « 2 • lO^^'^sec"^ when the object undergoing 
collapse is of size sa a. With such a small value of A, the "theoretical constraint" inequality (7.3) is violated: for 
the superposed sphere states considered in obtaining (7.3), the collapse time is w lOsec, much longer than human 
perception time. 

However, proponents of A = Ac could argue for a weaker "theoretical constraint" That is, when a human 

observer looks at the sphere, the detection process in the brain amounts to entangling the two spatially distinct sphere 
states with two spatially distinct states of brain particles. This extra entanglement, while only roughly estimable, 
appears to bring about collapse in less than human perception time. 

Detection of diffusion with such a small value of A could be possible. For example, we may compare the expected 
rotational diffusion (6.6) of standard quantum theory with the expected CSL diffusion (6.5) with A = Xq and a — 
IQ-^cm: AOqa/ « IQ-^irad and AGg ~ IQ-^i^/^rad. These give, for times (45min, 3hour), AOqa/ ~ (2.7, 10.8)rad 
and AGg ~ (1.4, 11.2)rad. 

Last, we call the reader's attention to Appendix F, where we consider that the collapse-inducing fluctuations of w 
may come from a thermal bath of some unspecified medium in thermal equilibrium with the 2.7°K cosmic radiation. 
The ~ t^/'^ time dependence of the CSL diffusion for a nucleon is identified with the standard Brownian motion at this 
temperature over an interval much less than the time it takes to reach thermal equilibrium. The latter time is taken 
to be 7-(the age of the universe), with 7 > 1. We then obtain the equality (F.2), A"^a^ w 10'^7, which is consistent 
with a wide range of parameter values, e.g., for a — lO^^cm this implies A^^ > 10^'^sec. 



VIII. EXPERIMENTAL CONSIDERATIONS 



In order to observe quantum mechanical rotational diffusion, either that arising from standard quantum theory or 
from CSL, it is necessary to isolate a small object from all outside torques for a period of minutes to hours while 
measuring its rotational position. We believe that it is now possible to perform such an experiment by combining 
techniques from nanomachining and from atom/particle trapping. Below, we shall consider the problems of creating 
suitable discs, suspending and isolating them, removing their residual thermal energy, and monitoring their angular 
position as a function of time. 

Our first consideration is the production of suitable disc samples. Larger discs (2/im diameter) have already 
been fabricated from silicon dioxide using standard IC fabrication methods and recent work at the Cornell 
Nanofabrication Facility shows that it is now possible to create structures with lateral dimensions below lOOnm and 
thicknesses less than 40nm [ p9[ . So it appears possible to make discs of suitable dimensions using current methods. 
As explained below, we suggest using discs of highly conducting metals such as copper or gold. 

We next consider suspending and isolating a disc for the duration of the experiment. We suggest utilizing a charged, 
conducting disc in a Paul trap in a very high vacuum. The Paul trap uses an alternating quadrupole electric field 
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to suspend a charged particle. While the method was initially developed to confine atoms, it was quickly adapted to 
suspend larger objects. Wuerker at. al. |^ injected small conducting microparticles using an electrostatic method 
that also charged the particles in the injection process. Once they had fed a cloud of particles into the trap they were 
able to select a single particle to retain in the trap by manipulating the trap's operating fields. 

More recently, Arnold and co-workers have used Paul traps and modified Paul traps to confine single micropar- 
ticles (a few /im in size) for optical experiments. One of their modified Paul traps has been used to confine single 
microparticles to within the Brownian limit set by the atmospheric gas in their traps. They add extra static electric 
fields to counter the effects of both gravity and imperfections in the quadrupole shape of the main field This 
leaves a perfectly force free spot in the trap where the particle will sit, subject only to collisions with the gas. Pressures 
of less than 5 • 10~^^ Torr have been reported in traps cooled with liquid He to 4°K. As we have remarked in 
subsection VIC, for a disc of radius 2 • 10^'^cm and thickness .5 • 10~^cm, at these conditions the average interval 
between gas-disc collisions is 45 minutes. Moreover, the Poisson statistical nature of the collisions make it likely to 
find intervals between collisions up to 90 minutes. This is quite long enough to observe even the rotational diffusion 
predicted by standard quantum mechanics. 

Although we are still engaged in studying the detailed dynamics of a charged conducting disc in the Paul trap, 
it appears already that the positional trap also acts as an orientational trap and will suspend the disc "vertically 
oriented" (with its flat surface parallel to the vertical symmetry axis of the field). There is no torque on a centrally 
positioned and vertically oriented disc causing it to rotate about the symmetry axis in the direction which we shall 
refer to as the azimuthal direction. Thus the trap appears to suspend a charged disc in exactly the best orientation to 
observe rotational diffusion in the azimuthal direction. What is currently under study is whether a displacement from 
center and/or tipping off vertical of the disc causes an azimuthal torque and, if so, whether that should be minimized 
or not (i.e., if the extra motion is due to translational or rotational CSL diffusion, this might cause an increase of 
observable CSL-induced diffusive rotation). 

When flrst injected into the trap, the discs will possess considerable translational and rotational kinetic energy. 
Arnold et. al. were able to remove this energy by the viscous interactions with the gas in the cell. In a high vacuum 
experiment there is no gas to take up this kinetic energy and a disc will continue to orbit the trap. However, if we 
add a transverse magnetic field then the eddy currents set up in the disc will convert the kinetic energy to thermal 
energy in the disc. The magnetic forces induced are proportional to the velocity of the disc and so provide a true 
viscous force. A simple dimensional analysis suggests that quite moderate fields, no more than a few kilogauss, will 
damp out the mechanical energy in a few seconds, thus bringing the disc to rest at the null point of the trap and 
vertically oriented. 

Preliminary calculations show that light (e.g., from a laser) shone in along the symmetry axis of the trap will scatter 
from the disc in a pattern that exhibits azimuthal anisotropy: more light is scattered perpendicular to the faces of 
the disc than perpendicular to the edges. Thus the orientation of the disc about a vertical axis can be monitored by 
collecting the scattered light. It appears that the geometry of the Paul trap makes it particularly easy to collect light 
scattered from a particle at the center of the trap. If the electrodes are highly polished, then the geometry is such 
that photons scattered away from the symmetry axis will be funneled by the electrodes to emerge through the two 
gaps where the cap electrodes and the ring electrode do not meet. Moreover, the scattered photons will retain their 
azimuthal orientation so that light collected at the gaps will retain the azimuthal intensity distribution and so provide 
information about the orientation of the disc. We suggest collecting the scattered photons with 8 photomultipliers 
operating as photon counters spaced around each gap so that the disc orientation can be measured to within 45°. 

The light which illuminates and scatters from the disc does so symmetrically and therefore exerts no average 
azimuthal torque on it. However, because the photons scatter randomly from the disc, they exert a random torque 
on it and so cause it to undergo diffusive rotation. Fortunately, this effect scales with the light intensity and thus can 
be minimized by using sufficiently weak illumination. Moreover, this illuminational diffusion can itself be measured 
in exactly the same way as any other rotational diffusion. Thus its effects can be eliminated by studying the behavior 
of the disc as the light level is reduced. The precise limit on the maximum amount of light that can be scattered 
without materially affecting the precision of the experiment depends on the value of A that one wishes to measure. 
For example, according to our present rough calculations, for the standard value of A = 10~^^sec~^, the disc can 
scatter about 200 photons per second before the illuminational diffusion exceeds 5% of the CSL diffusion. A simple 
calculation shows that, if the disc were at rest, you would need to count photons for about 10 seconds to localize it 
to within 45°. Since CSL diffusion with that A should lead to one revolution every 70 seconds the time resolution is 
quite adequate. 

In order to observe diffusion for a lower value of A, the maximum light level must be reduced accordingly. However, 
the lower value of A will lead to a slower rate of CSL diffusion and so allow integration of light over a longer period. 
This makes up for the lower maximum light level and means that the lower limit on A that can be measured is set 
by the vacuum and not by the illumination. A 45-90 minute interval between gas-disc collisions sets a lower limit 
A < 10^^'^sec^^. This is low enough to allow the experiment to definitively test CSL and, if the CSL diffusion does 



13 



not appear, to see the diffusion expected from standard quantum mechanics. 
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APPENDIX A: TRANSLATIONAL DIFFUSION IN CSL 

In CSL, the density matrix evolution of the wavefunction of a blob of matter containing N particles, in the position 
representation |xi, ...xat >= \x >, is given by 

8 A ^ ^ 

-(xIp(ON') - ~^{x\[H, PmW) - 2 ^ ^ + '^('^^ ^^-^ ^',)]{x\p{t)\x') (Al) 

i=i j=i 



where 



$(z) = e-^ , (A2) 



A is the collapse rate for a proton and we have assumed mass-proportionality of the collapse coupling (see section 
VII), Xj is the position coordinate of the jth particle, rrij is its mass, m is the mass of the proton and H is the usual 
Hamiltonian. In what follows we shall neglect the contribution of the electrons because of the smallness of the electron 
mass, and for simplicity take the mass of the neutron equal to m. Then, N is the number of nucleons. 

We wish to consider only the behavior of the center of mass (CM) of the blob. Accordingly, we trace Eq. (Al) 
over the relative coordinates Rj = — Q (eigenvalues r^) where Q = "^j-^*/ X]j "^j = is the CM 

position operator (eigenvalues q). $(xi — Xj) = $(ri — r^) and ^(x,^ ^ -"-j) ^ ^i^i ~ ^'j) ^^'^ independent of q but 
<I'(xi — Xj ) = <I>(ri — Tj- -|-q — q')- We shall also assume that the density matrix is the direct product of the internal and 
CM density matrices (this neglects their entanglement due to the collapse-induced excitation of the internal nuclear 
states). The trace of Eq. (Al) over relative coordinates yields 



d 

^(q|Pcm(t)|q') -j(q| 



p2 



N N 



A [ dr(r|p„t(i)|r)^^[$(r,-r,)-$(r,-r,-Hq-q')](q|Pcm(0|q')- (A3) 



= 1 3 = 1 



Since the nucleons are well-localized, we may write e.g., J dr{r\pint{t)\r)^{ri — Vj) « $(zi — Zj) where is the mean 
position of the ith nucleon. Moreover, since the nucleii are closely spaced compared to a = 10~^cm we may, to a 
good approximation, take them to be continuously distributed and replace the double sum in Eq. (A3) by a double 
integral, obtaining 



g^{q\pcmit)\q') = -i{q 



p2 



q') 

~^{v) I Iv ^^'^^'t*^^ " ^'^ - $(z - z' + q - q')](q|Pc™(t)|q') (A4) 

To see roughly how the collapse part of Eq. (A4) works, suppose that p(0) = (l/2)[|7/;i) -t- |V'2)][(V'i| + ('021] and 
that iV'i), IV'2) describe two well-separated (>> a) states of a blob so <i>(z — z' + q — q') sa 0. Therefore, neglecting 
the Hamiltonian term, Eq. (A4) says that the off-diagonal density elements exponentially decay: 



(q|Pcm(t)|q') = (l/2)(q|V'i)(^/^2|q')e 



-XNN't 



In this equation, if the dimensions of the blob are << a then <I>(z — z') « 1 and J dz = 1 so iV' « TV . If the 
dimensions of the blob are >> a, TV' « the number of nucleons in a volume a'^ . This collapse rate XNN' is diminished 
if the blob states overlap (<I>(z — z' + q — q') 7^ 0). 
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1. Translational Diffusion Of A Sphere 



We shall apply Eq. (A4) to an ensemble of spheres. Each sphere's CM wavefunction (subject to its own sample 
field w(x, i)) reaches an equilibrium size (sec Appendix B), subject as it is to the Schrodingcr evolution expansion 
and the collapse interaction contraction, with the center of a new contraction generally located off-center from the 
previous wavefunction center, thereby giving rise to the random walk. 

We shall use Eq. (A4) to calculate 



{Q^'){t) = / DwP{w) 



,{xi,,t\Qi'\i^,t), 



Tr[p{t)Q^ 



(A5) 



In Eq. (A5), \tp,t)w is the statevector of the sphere at time t which evolves under a specific collapse-causing random 
field w{x,t), the density matrix is p{t) = \4>,t)w yj{ip,t\ and, according to CSL, D'wP{w) = Dwyj{tjj,t\'il),t)yj is the 
probability that the field w{x, t) appears in nature, where Dw ~ Hx t dw{x, t) (space-time may be regarded as divided 
into little cells, in each of which t«(x, t) can take on any real value). In Appendix B we shall calculate 



{Qifit) = J DwP{ 



w) 



,{^P,t\Q^\ij,t}, 



(A6) 



which cannot be expressed as a trace with respect to the density mat rix. As shown in Appendix B, the spheres we 
consider arc large enough so that the mean square packet width = {[Q^ — (Q^)]^) rapidly reaches an equilibrium 

constant size << {Q^)^{t). Therefore the increase with time ~ of {Q^)^{t) found here (see Eq. (AlO) is solely due 
to the diffusion of the spheres. 

To find {Qi'^){t) we take successive traces of Eq. (A4): 



(A7a) 



(A7b) 



d 2 XN^h'f{R/a) 



(A7c) 



where 



f{R/a) = 



1 



dzrfz'$(z - z') 



1 - 



2a? 



Integration of / may be facilitated using Gauss's law to convert the volume integrals to surface integrals: 
/(i?/a) = ^ j J dzdz'V -ejV' -e'j^iz-z') =2a^^ J J dA ■ ejdA' ■ e'j^{z - z') 



(A8) 



(A9a) 



2a^ 



2a^ 



(A9b) 



R« a 



/(I) = -62, 



R» a \R 



(A9c) 



/ is a monotonically decreasing function of its argument. 
It foUows from Eqs. (A.7) that 



{Qn = {[Q' + 



PH 



)(0) + 



Xh'^f{R/a)t^ 



(AlO) 



15 



which is the result quoted in Eq. (3.4). The difTusion term in Eq. (AlO) can be understood as follows. cP{Q^^)(t)/dt^ 
is proportional to the square of the collapse-induced velocity {h/Ma)'^ multiplied by the effective collapse rate. For 
R << a, the collapse rate is ^ XN'^, giving rise to Eq. (AlO) with f — 1. For R >> a, as we have previously 
shown 1^, the collapse rate is ~ A- (number of particles in a volume a^) ■ (number of uncovered particles). That is, we 
imagine the sphere in a superposition of two states displaced from each other by a certain distance so the two images 
of the sphere overlap: the "uncovered" particles are those in the region of no overlap. In this case we suppose the 
displacement distance is a. Then, the (number of uncovered particles)^ {N/V)a{the surface area A of the sphere). 
Thus we find the expression 

- {h/MafX{Na^/V){NAa/V) - \{n / maf {a / Rf - X{h/maff. 



2. Translational Diffusion Of A Disc 



In the case of a disc undergoing translational diffusion, / depends upon its orientation. If the disc is of radius L and 
thickness b, for motion perpendicular to the disc face it follows from Eq. (A9a) (which is applicable to an arbitrarily 
shaped object) that 



/ = 4 

For example, for (6/2a)^ 
the disc edge, 

/ 



r, \ 4 /„ \ 2 

2a \ I 2a 



L/2a 



xdx 



L/2a 



(All) 



<< 1, / w 1 for (L/2a)2 << 1 and / {2a/Lf for {L/2af » 1. For motion parallel to 



- e /,(L/2a) - 



b/2a 

— I dxe~ 



b 

2a 



-{hl2af 



b/2a 



(A12) 



where / ^ (4/v^)(a/L)3 for {b/2af « 1 and {L/2a)^ » 1. 

Eqs. (All, A12) can also be understood as proportional to the effective collapse rate. We use an alternative 
expression for the rate, equivalent to that given in the previous paragraph, appropriate for an object with a dimension 
{b in this case) less than a. It is rate ^ (number of particles in a cell) ^-(number of uncovered cells), where a cell is a cube 
of dimension a on each side. Here each cell has occupied volume ba^ so the number of particles/cell = {Nba^ / nL'^b) ^ 
[a/L)"^. For displacement a » b perpendicular to the face, all the cells — i^l? ja^ of them — are uncovered, so we 
obtain the collapse rate ^ {a/ LY{L / a^ = {a/L)'^. For motion parallel to the face, displacement a uncovers the cells 
lying on the circumference of the disc ~ 2TrL/a of them, giving the collapse rate ~ [a/ L)'^{L/a) = (a/L)^. 



APPENDIX B: WAVEPACKET WIDTH OF CENTER OF MASS IN CSL 



The CSL evolution equation for the normalized statevector in Stratonovitch form (so manipulations can be per- 
formed using the usual rules of calculus) is 



= <^ -iH 



dxG(x)w(x,t) - A J dx[G^{x)- ^{^|J,t\G^{x)\tP,t), 



(Bla) 



G(x) 



AT p 



2^3/4 



'(Xj-x)^/2a^ 



,{ip,t\e 



^)/2a^ 



(Bib) 



where Xj is the position operator for the jth nucleon, w{Ti,t) = dB{x.,t)/dt is standard white noise and B(x,t) is 
standard Brownian motion {w{x, t) = 0, w(x, t)w{x' , t') — A(5(x — x')6{t — t')). We extract the equation for the CM 
wavefunction just as in Appendix A whose notation is used here (Eq. (Al) for the density matrix can readily be 
derived from Eq. (Bl)). Again, we assume that the statevector can be written as a direct product of the internal 
statevector \tpint,t) and the CM statevector \(j),t)i^. We suppose that |'0int,i) obeys the usual Schrodinger equation 
(thereby neglecting the CSL excitation of atoms and nucleii) with Hamiltonian Hint , so the complete Hamiltonian is 
H = P'^ /2M + Hint, where P is the CM momentum operator. Using this in Eq. (B.l) with Xj = R_, +Q, multiplying 
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by J dr{il;int,t\r){r\, employing the localized nature of nucleons so, e.g., J dr\{il;int,t\r)\'^ F{rj) « F{7,j) where Zj is 
the mean position of the jth nucleon, and then approximating F{zj) « (N/V) Jy dzF{z) results in 



dt 



(q|'/',i)t 



|-«^+ Jdxg{x-ci)w{x,t)-xjdx[g'^{x-ci)-,j„{i;,t\g'^{x-ci)\-tp,t)y,] (B2a 



^(x - q) 



dz 



,t\e- 



(7ra2)3/4 \V J Jy 

(Eq. (A4) for the CM density matrix can readily be derived from Eq. (B2)). 



(B2b) 



1. A Sphere's Equilibrium CM Wavepacket Width And The Time To Reach It 

Eq. (B2), applied to a sphere of radius R, is the starting point for our calculation. We shall consider only cases where, 

for each \(p,t)w, the squared wavepacket width s'^{t) = w{(l),t\Q^'^\(j),t)w— w{<l},t\Q''\4'jt)-ui'^ = w{(f'^t\[Q'' — {Q'')Y\<l>^'t)w 
is much less than (note that s has no subscript j because we assume its initial spherical symmetry which is 
maintained thereafter). At the end of section Illit is shown that s « a implies R > 10~^cm. 

We expand the exponents in Eq. (B2b) in powers of [Q^ — {Q^)]/a, retaining only the leading term: 



- - (y) X dze-^'^<^^-r/^'^-}^, + (Q) - X) . (q - (Q))/a^. 

The solution of Eq. (B2a), for a short time At, can be written (with use of Eq. (B3)) as 



(B3) 



{q\4>,t + At)u, = exp 



N' 



-XAt^f{R/a)[{ci-{Q)f - ^(V,t|(q- (Q))'|^,^)^ 



{fi.\<i>,t)z 



(B4) 



(note the replacement of x — (Q) by x as dummy integration variable and the concommitant use of translation 
invariance of w{x,t)) where f{R/a) is given by Eq. (A8). Eq. (B4) shows that a gaussian wavefunction at time t is 
taken into a gaussian wavefunction at time t + At. 

Although we could deal with a more general class of wavefunction, the results are the same and the argument is 
simpler if we restrict ourselves to the complex gaussian wavefunction 



(B5) 



In Eq. (B5), b = b/{ + ib/, a'^ = a\+ iaj are complex functions of time. Using this wavefunction one can calculate 
various expectation values involving the CM position and momentum. It follows from Eq. (B5) that 



|(q|</.,t)j2 = (27r.2)-V4^-(q-(Q>W, 
{Q)=hn + hjaj/a%, (P) = b,24, ^ {{Q^ - {Q^m = aj, + /aj,, {{P^ - {P=)f) = l/Aal- 

We note that 

Putting this into Eq. (B4) and equating the coefficients of q^ and q results in 

1 d 1 i XN"^ 



(B6) 



4dta^ 



SMo-^ 2a2 



■/ 



(B7a) 



l_d b_ 



ih 



First, consider Eq. (B7a): 



N 



a2(7ra2)3/4 y 



J dz J dx(z 



x)e' 



-(z-x)^/2a' 



w{x, t) 



(B7b) 
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2XN^ 



-fa' 



(B8) 



It has no stochastic part and may be immediately solved. As t ^ oo, where da'^{t)/dt = 0, according to Eq. (B8), 

£7^(00) = (a/2Ar)(n/2MA/)i/2(i + i) 
so, by Eq. (B6), the asymptotic squared wavepacket width is 

s2(oo) = sl, = {a/N){n/2M\fY/\ (B9) 



This result can be obtained by a simple physical argument given in section III (following Eq. (3.5)). 
Introducing 0-^(oo) into Eq. (B8), together with 



converts Eq. (B8) to 



(BIO) 



with solution 



d{tlT,) 



a 



£!W [e*(i+^)/-» + 1] + (l+l)[e*(i+0/r, _ 1] 

£fW [gt(l+i)/r3 _ 1] + (i+i) [et(l+i)/r. + 1] 



(Bll) 



(B12) 



which shows the approach to equilibrium. 

Thus we have achieved the main purpose of this appendix, to obtain the cm wavepacket equilibrium width (B9) 
and the characteristic time to reach that width (BIO). We emphasize that these results apply to every cm wavepacket, 
since they are independent of the particular realization of the fluctuating field w encountered by a sphere. 



2. CM Translational Diffusion Revisited 



However, we do have Eq. (B7b) which does depend upon w and which gives us, via Eq. (B6), each individual 
cm wavefunction's mean position (Q) (and mean momentum (P)), enabling us to understand in detail the ensemble 

average {Qi^} in Eq. (AlO). 

We first note that the stochastic term in Eq. (B7b) has no dependence on the dynamical variables ct^ and b. Since 

the ensemble average of the product of this term's ith and jth components is SijXN"^ fS{t — t')/2a^ we may write the 

term as (AA^^//2a^)^/^w(t) where the Wj(tys are independent white noise, 'Wi(t)wj{t') = 5ijS{t — t'). We may then 
write Eq. (B7b) as 

4b = -^^ + ^. (B13) 

Suppose we follow an individual wavefunction for sufficient time >> until (say, at time t = 0) it achieves its 
equilibrium width Sqo with cr^ = s^(l + i)/2 so Eq. (B13) simplifies to 

dh = ^dt + il±i) ^dB(t) (B14) 

where B{t) is Brownian motion (w(t)=dB(t)/dt). The solution of Eq. (B.14) is 

= i^J' + i^^BKt), m = i^.B^it) (B15) 

(we have assumed 6jj(0) = bj{0) = 0). 
It follows from Eqs. (B15) and (B6) that 
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(Q) 



o 3/2 , 

2ts' Jo 



Jo Tc 2s^r«' 



(B16) 



which expHcitly shows the diffusive nature of (Q) and (P) 



We can now find (Q-'^) 



(Q^)^ and compare with Eq. (AlO). RecaUing that B:i\t) = t and 



B3{t)B3{t') =min(i,f'), we obtain 



t 

+ 



+ 



which is identical to Eq. (AlO) for a wavefunction which has equihbrium width at i = 0. 



(B17) 



APPENDIX C: ROTATIONAL DIFFUSION IN CSL 



Starting with Eq. (Al) for the evolution of the density matrix in CSL, we follow the lines of argument in Appendix 
A. We assume here that the cm of a blob of matter is fixed but that it is free to rotate about a fixed axis through an 
angle represented by the operator O with angular momentum operator C. We also assume that the density matrix is 
the direct product of the internal density matrix and the orientation density matrix Pang- We obtain, analogous to 
Eq. (A4), 



d_ 
di 



{9\Pang{t)\e') = -m 
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) Pang {t) 



\e') 



-A 



J J^dzd^'mzio) - z'(o)) - $(z(0) - z'{e'me\pANG{t)\e') 



where, denoting the rotation axis by Z3, 



1 



^{z{0) -z'{e')) =exp-^[z^ +z'2 - 2{ziz[ + Z2Z2) cos{9 - e') - 2(zi4 - Z2z[) sm{9 - 9') - 2zzz'^. 



To find (©^)(t), in analogy to Eqs. (A7), we take successive traces of Eq. (CI): 

d 



(Cl) 



(C2) 



(C3a) 



|i(/:e + e/:) = ^(/:^) 



(C3b) 



--W) = - 

dtP^ ' 2 



Irot 



where z±_ = {zi, Z2) and fnor is the dimensionless geometrical factor 



(C3c) 



fnoT = 2 



Ma 

w 



dzdz'[z^ ■ z'_L - ^(z-L X z'_l)2]$(z- z'). 



(C4) 



To see that (C3) vanishes if the blob is rotationally symmetric about the 2;-axis (i.e., a sphere or a disc with the 
z-axis perpendicular to its face), we write (C4) as 



fnoT 



IV 



> 

[ dz{z^xV^J^ [ (iz'$(z-z'). 
Jv Jv 



(C5) 



Rotational symmetry implies that the integral over z' is just a function of z^ and then the integral vanishes since 

(z_L X Vz_l)z^ = 0. 

It follows from Eqs. (C3) that the mean square angular diffusion has the time dependence 
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(e^) = ((e + y)')(o) + A^ 



1 2 



ROT- 



(C6) 



We wish apply Eq. (C6) to a disc of radius L and thickness b {I = (AfL^/4)[l + (&^/3L^]), with the rotation axis 
parallel to the face of the disc, so Jrot rnust be calculated for this case. The double volume integral in Eq. (C5) can 
be converted to a double integral over the surface of the disc by using the divergence theorem: 



/i?,OT(a,/3) 



[l + (/32/3a2)]/3a4 



[ dA • (r X k) [ dA' ■ (r' x k)e^ 

J A J A 



ir-r'f 



(C7) 



where k is the unit vector along the axis of rotation and a — L/2a, /3 = b/2a. Calling the contribution of the two 
disc faces /i, the two disc edges /2, and the edge- face contribution /a, we obtain 



/rot {a, (3) 



[1 + (/32/3a2)]/3a4 



[/i + /2 + /s] 



(C8a) 



Jo Jo 



(C8b) 



h = (l/2)a2e-2" /^(2a2) / / y'dy'e-^y-y > 

J-/3/2 J-0/2 



(C8c) 



ra r0/2 

/3 = -2ae-"'/ r'^dre^''^h{2ar) ydye~^-'^'^-y^' 

Jo J-P/2 



(C8d) 



where Ji is the Bessel function. 

A graph of fRorioi, (3) vs a parametrized by various values of [3 is given in FIG. 1. We note that, for a thin disc 
(/3 << a) with (3 « 1, /i is the leading term in Eq. (C8a) which becomes 



/rot {a) 



r^dr / r'2dr7i(2rr')e-('-'+'^"). 

-'0 



(C9) 



APPENDIX D: THERMAL RADIATION VISCOSITY FACTOR FOR A DIELECTRIC SPHERE 

In order to compare the Brownian diffusion of an object in a thermal radiation bath with CSL diffusion, it is only 
necessary to find the viscosity factor ^ for this situation and put it into the Brownian motion equations of section 
II. An object moving with respect to thermal radiation with speed v feels a drag force — because it receives more 
momentum from the photons it approaches than from those from which it recedes. We have not been able to find ^ 
for a dielectric sphere in the literature (the closest has been the force on an oscillator Q|) so we give it here. Actually, 
after this Appendix was written, we decided that the experiment we propose would concern a conducting disc rather 
than a dielectric sphere! The result for a conducting sphere is not quite the same as that for a dielectric sphere with 
dielectric constant equal to infinity: although that is the appropriate limit for electric field behavior, a conducting 
sphere's magnetic behavior is also important in considering the scattering cross-section of electromagnetic radiation 
(necessary for this calculation). And, of course, a sphere is not a disc. However, the result obtained for ^ will be 
representative, i.e., the same up to a numerical factor not too far from 1, when the dielectric constant goes to infinity 
and the radius of the sphere is replaced by the radius of the disc. 



1. Viscosity Factor For a Mirror 



For expositional ease and purposes of comparison we shall first obtain Einstein's result for a mirror moving perpen- 
dicular to its face |21^J23], as seen from the laboratory frame in which the radiation is thermal. We shall use properties 
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of photons (which Einstein had not yet obtained as he was in the process of estabhshing these, so he used classical 

electromagnctism) . 

As is well known, at temperature T the mean number of photons in a mode of frequency (the subscript refers 
to the laboratory frame) is [exp f3hvo — (where (3 = {kT)~^). Since the number of photon modes/vol of frequency 
vq in the range dvo moving in a direction (^O) </'o) within solid angle dflo is 2pQdi'odQo/c'^ (the factor 2 is for the two 
polarizations), the mean photon number /vol- freq-solid angle is 

n(i/o) = 2(i/2/c3)[exp/J/iz/o-l]-' (Dl) 

First we find the momentum transferred to the mirror by a colliding photon. Einstein considered a mirror which, 
in its rest frame, is perfectly reflecting only for frequencies in the range {v, v + dv) (no subscript refers to the rest 
frame of the mirror) and perfectly transmitting otherwise. Let the mirror (of area A) move in the ^-direction with 
speed V away from a photon of momentum po = hv^ / c whose direction of motion makes an angle 9q with respect to 
the 2-axis. The photon will be reflected only if it has frequency v in the rest frame of the mirror. From the energy 
and momentum transformations of special relativity (all calculations are to order v /c) , 

V = vq[1— [v/c)cos9o], vcosO = vqIcosOq — [v/c)] (D2) 

where 6 is the angle the photon makes with the normal to the mirror in the mirror's rest frame. In this frame the 
photon's incident and outgoing (normal) momenta arc respectively {hv/c) cos 9 and —(hv/c) cos 9 so the momentum 
imparted (normal) to the mirror is AP = 2{hv / c) cos9 . The difference of momenta of a nonrelativistic object is a 
Galilean invariant. Therefore, AP = APq which, by (D2) may be written as 

APo = 2{hvQ/c)[cos9o-{v/c)]. (D3) 

Next we find the number of these photons colliding with the mirror in time dt. In the mirror rest frame this is 
J • Adt where J is the particle number flux and A = At, with A the area of the mirror. This is the same number 
which collides with the mirror in the laboratory frame in time dt. The four-current transformation equation gives 

J • A = (Jo — pov) • A = n{vQ)dv(id,Vio{c cos 9^) - v)A. (D4) 

Thus, by Eqs. (D3, D4), the momentum transferred to the mirror in the laboratory frame in time dt, expressed in 
laboratory frame coordinates, is 

- vdtdS, EE J • AdtAPo = n{vo)dvodVlnc[cos9n - {v j c)\Adi2(hvQl c)\cos9q -~ (u/c)]. (D5) 

It remains to integrate Eq. (D5) over all fio but, first, we must express vq in terms of v and 6'o- From the inverse 
of Eq. (D2) we have vq = v\\.-\- {v/c) cos^o] so we obtain 

dvQ = dv\l + iv/c) cos6'o], voniyo) = H — cos9Q]r]\ u\l -\ — cos9q\ ) = vniv) + (vniv))' viv / c) cos9q + oiv/cf' . 

c \ c J 

Then, we must remember that the above analysis is predicated upon the mirror receding from these photons (so the 
range of 9o is (0, 7r/2)). The momentum imparted by the photons on the other side of the mirror is given by the 
negative of the right hand side of Eq. (D5) with the replacement v — > —v. Thus, the contribution from all photons 
to the force is 

- vd4= dv2hA j rfOo|z^«(z^)[cos^ - (w/c)(2cos6'o - cos^ ^o)] + {vn{v))' v{v / c) cos^ 
— {v — > —v) 

= -vdv2-K{h/c)A[2,un{u) - u{un{u))']. (D6) 
This is Einstein's result. Putting Eq. (Dl) for n(i') {v = uq to zeroth order in v/c) into Eq. (D6) yields 

Of course, v may be integrated over to obtain the viscosity factor for a mirror which is a perfect reflector at all 
frequencies: 
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2. Viscosity Factor For A Dielectric Sphere 



Our discussion for a dielectric sphere (dielectric constant e, radius R, moving in the z-direction with speed v) 
exactly parallels that for the mirror. 

First we find the momentum transferred to the sphere by a colliding photon. In the rest frame of the sphere, the 
scattered radiation has a dipole pattern so the radiation scattered in two opposite directions carries no net momentum. 
Thus, for radiation of frequency j/, insofar as momentum transfer is concerned, the sphere acts like an absorber (of 
area equal to the total scattering cross-section cr(i^)). Thus, the momentum effectively imparted (i.e., on average) in 
the z-direction by an incident colliding photon is AP = {hv/c) cos9. As in our previous discussion, since AP — APq, 
the effective momentum imparted by a single photon in the laboratory frame is 1/2 of the value given in Eq. (D3). 

Next we find the number of these photons colliding with the sphere in time dt. In the sphere rest frame this is 
J(7{i')dt, where J is the number flux: this is the same number that collides with the sphere in the laboratory frame 
in time dt. To express this number in terms of laboratory frame variables, we note that J/c is the zeroth component 
of the current 4- vector and J cos 9 is the component along the z-axis. Therefore the Lorentz transformation of the 
zeroth component of the current 4- vector is J/c = [Jq/c — {v/c^)Jq cos^q] or, substituting for Jq, 



J(j{i')dt ~ n{vQ)dvQdVLQc[l — (v/c) cos6Q]a{i')dt. 



(D9) 



Note that Eq. (D.9) differs from the parallel mirror equation (D.4) in that radiation of frequency incident from any 
direction sees the same cross-section of the sphere while this is not the case with the mirror. 

Therefore, the momentum transferred in the z-direction in the laboratory frame in time dt by these photons is, by 
(D9) and half of (D3), 

— vdtd^ — Ja{v)dtAPQ — n{vo)dvod^lQc[l ~ [v/c) cos 9Q]a{v)dt{hvf) / c)[cos9q — {v/c)]. (DIO) 
As before, we express dv^ and n(i/o)j^o in terms of v and 9q {a is already in terms of v) and integrate over all Qq: 



vdS,— dvha{v) i dQ,Qlvn{v)[cos6(i — {v/c)\'V{vn{v))'v{v/c)cos^9Q 



-vdvi-iTT /'i){h/ c)a{v)[ivn{v) - v{vn{v))'] 



(Dll) 



This is 2/3 of the comparable expression (D6) for the mirror, with the area A replaced by the cross-section cr{v). 

The classically calculated cross-section (i.e., the total scattered energy/sec divided by the incident energy/sec-area) 
for an electromagnetic wave of wavelength > > i? is W5[ 



87r\ f2'Kv 



.WHyji — 



(D12) 



where, for simplicity, we shall only use the limit of large e. In Eq. (D12), a has been averaged over incident 
polarizations and summed over scattered polarizations. 

Putting (Dl) for n{v) and (D12) for a{v) into (Dll) yields 



f hi 



dv. 



(D13) 



We remark that, if dv(j{i') in the first equation of (D13) is replaced by Jp°° di'a{y) = ttc^ /mc^ the sum rule for an 
individual oscillator | |4^ of mass m and resonant frequency v, we obtain the value of ^ for a single oscillator given in 
reference [Q. 

Upon integrating D13 over v we obtain the viscosity coefficient 



dzz^ 



4(2.)V^V,^a 



[e^ - 1]2 135 \hc J 



(D14) 



since the integral=(27r)^/60 (sa 8!). This result is used in Sections IIB and IVC. 
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APPENDIX E: A GRAVITATIONAL PROPOSAL 



Diosi Ql suggested a gravitationally based CSL-type collapse model with the collapse rate ^ G. However, it 
effectively had a w the proton size and therefore too large a proton excitation rate, a flaw corrected by Ghirardi, 
Grassi and Rimini who added the standard a to the model. 

Penrose Q, perhaps unwilling to commit to a nonfundamental parameter a (however, see Pearle and Squires [3^ ] 
for a "derivation" of a in terms of fundamental constants in the context of a gravitationally based model) has a 
more modest proposal. His suggestion is that, when quantum theory describes an object as being in a state of two 
superposed locations, collapse of the state to one of those locations will take place in a time equal to fi divided by 
the gravitational energy required to move two real copies of the object from a completely overlapping configuration 
to these two locations. For example, consider a sphere of mass M and radius R. Since the gravitational energy of two 
such spheres displaced by a small distance D << R is 

U{D) = 

then the time it takes a quantum state of a sphere in a superposition of two states separated by the distance D to 
collapse to one or the other state is 

Tc = 2hR^/GM^D^ (El) 

This is a minimalist proposal, not a complete dynamical theory. For example, it is silent on how to treat the collapse 
of the state of a sphere in a continuous superposition of locations (i.e., the usual wavefunction description of the CM 
of a sphere). Nonetheless, we shall have the temerity to make what we regard as a reasonable extrapolation to that 
situation, in order to estimate the random walk entailed by this proposal. 



GAP 


6 






R 


^ 5 ^ 







1. Equilibrium CM Wavepacket Size For A Sphere 

First, consider the qualitative argument given after Eq. (3.5), for the equilibrium size of a CM wavefunction, applied 
to the sphere. A CM wavepacket of width D expands a distance ~ {T%/MD)lS.t in time Ai due to the Schrodinger 
evolution. Now, assume that the collapse is linear, in the sense that, in time At, if the wavepacket width is D, collapse 
acting alone makes it contract to D[l — {/Si / t^)] , where is given by Eq. (E.l). If D = s is the equilibrium width of the 
wavepacket, then the Schrodinger expansion is compensated by the collapse contraction, yielding sAt/rc {h/Ms)At 
or 

' GM^ 

Eq. (E.2) may be compared to the CSL result (3.5): 

s - ' 



XM3f{R/a)' 

We may therefore regard this proposal's equilibrium CM wavepacket size as giving the CSL size if 



Gm^ f a^^ 



In particular, if i? ~ a (and so / ~ 1), 



A ^ ^ « m-^hec-' (E4) 
an 
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2. Translational Diffusion Of A Sphere 



We may obtain the same result, that this proposal gives the CSL behavior for objects of size « a with A having 
numerical value (E.4), from other considerations such as random walk of the sphere. In this case, the Schrodinger 

equation tells us that d^{Qi'^)/dt^ = {2/M'^)d{Pi'^) /dt (Eqs. (A7)). The collapse, acting on a wavefunction of width 
D, narrows the wavefunction and, in so doing, increases the energy. As we have remarked in the previous subsection, 

dD/dt = -D/tc = -GM^D^/2hR^. 

Prom the uncertainty principle, (P-'^) ~ {h/D)'^, so 

d{P^^)/dt ~ -{h^/D^)dD/dt ~ GM'^n/R^ 
(notice that the result is independent of D, as in CSL) and so 

d^'iq^/dt^ ~ Gn/R^ (E5) 
(notice that the result is independent of M, as in CSL). Comparison of Eq. (E.5) with the CSL result (3.4): 

d^W)/dt' = ^^"fWa) 

yields the same "effective" A given in (E.3) 



3. Rotational Diffusion Of A Disc 



Angular random walk of a disc proceeds along the same lines. For a thin disc of mass M and radius L, the 
gravitational energy required to rotate one such disc through a small angle 9 with respect to a second initially 
completely overlapping disc is ~ {GM"^ /L)6'^ so 

Tc-^hL/GM'^e'^. (E7) 

Here we utilize, from Eqs. (C.3), d^{9'^)/dt^ = {2/P)d{£?)/dt. According to this gravitational proposal, d6/dt = 
—0/tc- From the uncertainty principle, {C^) <^ so 

d{c^/dt ~ -n^/e^de/dt ~ n^/e'^Tc 

and so 

d^le^/dt^ ~ h^/I^0\e ~ Gh/L^ (E8) 
Eq. (E.8) may be compared with the CSL result (6.5) 

d^Je^/dt^ ~ \{n/ma'^ffROT{L/2a) (E9) 

which yields the "effective" A 

AW(i/2a)~^(^|j (ElO) 

In our proposed experiment, for which L a (and so fiiOT{L/2a) « 1), the "effective" A is again given by (E.4). 
The results obtained here are effectively the same as would be obtained with the modified Diosi model. 
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APPENDIX F: THERMAL SOURCE OF THE FLUCTUATIONS? 



It is fun to speculate that the collapse-inducing fluctuations of w may come from a thermal bath, as do so many 
other fluctuations in physics. 

Since a thermal bath defines a preferred reference frame (i.e., the frame in which the bath medium has zero average 
momentum density), this would preclude a truly special relativistically invariant collapse model. But, anyway, the 
universe is not truly special relativistically invariant, possessing as it does the preferred comoving reference frame. 
Moreover, this reference frame is endowed with the 2.7°K thermal radiation bath. So, one might entertain the idea 
that the fluctuations of w arise from some unspecifled medium in thermal equilibrium with the 2.7°K radiation. 

For an object in random walk, we note that the ^ time dependence of (Q-'^) given by Eq. (3.4) for CSL is also 
the time dependence of (Aa;)^ given by Eq. (2.7c) for ordinary Brownian motion when t is very much smaller than 
r = ^/M (which characterizes the time scale of the approach to thermal equilibrium). Since objects show no sign of 
reaching thermal equilibrium today, we may assume that r is larger than the age of the universe, r = 750A^^j^ (since 
A^^^ = IQiSecw 3 • 10*yr) with 7 > 1. 

Continuing in the same lighthearted vein, we propose that, for a fundamental object, the nucleon, the two sources 
of the ~ behavior, thermal and CSL, may be identified, and we equate Eqs. (2.7a) and (3.4), obtaining 

« „ or KT « 5O7 „ „ (Fl) 

(we have set M = m and f{R/a) = 1). 

When T = 2.7°K then kT « 2.5 • IQ-^eV. The energy h^/m?a%s^ « 4 • IQ-^eV. Thus, Eq. (F.l) implies 7 ^ 10^, 
which is consistent. 

Of course, in the speculation above there is no need to choose (A~^,a) to have their CSL numerical values. The 
appropriate generalization of (8.16) is 



A ( a \ 7 



-j-r or A-^a^ w 10^7. (F2) 



\Q\j^\acsLj 10^ 

There is quite a range of A and a consistent with (F.2) and present constraints (see FIG. 2), especially in view of 
the flexibility in choosing 7. 
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FIG. 1. A graph of fROT{a,f3) as a function of a s(disc radius)/2a, for various values of (3 s(disc thickness)/2a. 



FIG. 2. A graph of logj^g A ^ vs logj^g a for various constraints: the constraint boundaries specified in Eqs. (8.1), (8.2), (8.3), 
(8.4) are respectively labelled 1, 2, 3, 4. 



TABLE I. CSL diffusion in vacuum: rms distance AQcm for various sphere radii R and times t. 



t in sec 



R in cm 




10 




10^ 




10^ 


lO"" 


8 


10"'^ 


8 


• 10"^ 




8 


10"' 


6 


10"" 


6 


• 10"^ 




6 


10-'' 


2 


•10-^ 


2 


• 10"* 


2 


•10-1 


10"^ 


6- 


10-" 


2 


• 10"^ 


2 


• 10-^ 


1 


6- 


10-15 


2- 


10-12 


2 


•10-9 
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TABLE II. CSL rms equilibrium center of mass wavefunction size Soo and characteristic time Ts to reach that size in vacuum 
for various radii R oi a sphere of density Igm/cc. 



R in cm 


Soo in cm 


10"'' 


7- lO"'' 


10"^ 


4 • lO-'^ 




1 • 10"® 


10-2 


4-10"" 


1 


1 • IQ-" 
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